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The spontaneous magnetization as a function of temperature is investigated for a number of 
disordered Ni-Fe-Mo and Ni-Fe-W alloys using superconducting quantum interference device mag- 
netometry, with a focus on the low-T behavior as well as the critical exponents associated with 
the magnetic phase transition. While the low-T magnetization is found to be well described by 
Blochs T 3//2 law, an extraordinary enhancement of the spin-wave parameter B and the reduced 
coefficient B 3 / 2 =BT^ 2 are observed with increasing Fe dilution as compared to conventional 3d 
ferromagnets, whereas the critical amplitudes are found to decrease systematically. Recent locally 
self-consistent calculations of finite-temperature spin dynamics in a generic diluted magnet provide 
an understanding in terms of two distinct energy scales associated with weakly coupled bulk spins 
in the FM matrix and strongly coupled cluster spins. In view of similar behaviour observed in di- 
luted magnetic semiconductors and other ferromagnetic alloys, it is proposed that these distinctive 
features corresponding to the three important temperature regimes provide macroscopic indicators 
of signature effects of spin clustering on magnetization behaviour in disordered ferromagnets. 

PACS numbers: 75.30.Ds, 75.40.Cx 



I. INTRODUCTION 



Magnetic phases of alloy systems with both ferromag- 
netic and antiferromagnetic exchange interactions have 
been studied in details by many groups including ours. 
Ni-rich Ni-Fe-Cr and Ni-Fe-V alloys^ 2 - have shown ex- 
otic magnetic phases and hence have been probed ex- 
tensively. Compositions near Ni80% - Fe20% have been 
widely investigated in bulk form not only for their rich 
magnetic phases but also for the unique quality of low 
magnetostriction and high initial permeabilitj* 3 -. These 
are ideal properties for making magnetic cores for light 
electrical equipment as well as magnetic shielding mate- 
rials. 

The observation of several signatures of a spin-glass- 
likc phase in recent dc magnetization and ac susceptibil- 
ity studies of Ni-Fe-Mo and Ni-Fe-W alloys^— indicates 
presence of competing ferromagnetic and antiferromag- 
netic spin interactions, resulting in frustration-induced 
frozen spin configurations in spin clusters. Generally, a 
broad distribution of magnetic spin couplings would be 
expected corresponding to different separations and con- 
figurations resulting from the positional disorder of mag- 
netically active Fe atoms in these alloys. 

Are there any distinctive finite-temperature spin dy- 
namics effects associated with this disorder-induced dis- 
tribution of magnetic interactions in the Ni-Fe-Mo and 
Ni-Fe-W alloys which are observable in their macro- 



scopic magnetization behaviour? A similar question 
was recently addressed in the context of squid mag- 
netization studies^ of diluted magnetic semiconductors 
Gai-a-Mn^As, which also exhibit strong disorder effects 
due to Mn positional disorder. Recent theoretical studies 
of finite-temperature local spin dynamics within a generic 
model for diluted magnets have provided fundamental 
understanding of macroscopic magnetization characteris- 
tics in terms of microscopic spin disorder, spin clustering, 
and distribution of spin couplings^— 

There are three important temperature regimes in the 
magnetization behaviour of a generic disordered ferro- 
magnet. The low-temperature regime characterized by 
the spin- wave parameter, the intermediate-temperature 
regime nearly upto Tq characterized by the overall shape 
of the magnetization decay, and the critical regime very 
close to Tq involving critical fluctuations, divergent spin 
correlation length, and characterized by the critical expo- 
nents and critical amplitudes. In this paper we propose 
that the magnetization behaviours in all three tempera- 
ture regimes actually provide macroscopic indicators of 
signature effects associated with microscopic spin dynam- 
ics of weakly coupled bulk spins and strongly coupled 
cluster spins. For that purpose we have made a detailed 
study of the temperature dependence of magnetization 
in Ni-Fe-Mo and Ni-Fe-W alloys and the critical expo- 
nents and critical amplitudes associated with the phase 
transition. 
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II. EXPERIMENTAL DETAILS 

Alloys of both ternary series (Mo and W) were pre- 
pared by arc melting of the required amount of con- 
stituents of Johnson-Mathey, 99.999 % pure Ni, Fc, Mo, 
and W in pure argon atmosphere. Then the shiny but- 
tons were sealed in quartz tubes, flushed with ultra pure 
argon gas, and homogenized at 1300°C for 48 hours. 
Then they were cold-rolled and cut to size. A final an- 
nealing was done in argon atmosphere at 1100°C to re- 
duce strain due to cold rolling. We kept the concentration 
of Ni roughly around 80 - 83 at. % and varied Fe from 17 
to 3 at. % and Mo from 2 to 14 at. %. For alloys contain- 
ing W, Ni concentration remains roughly around 79 to 86 
at. % while that of Fe steadily decreases from 18 to 3 at. 
% and that of W increases from 3 to 10 at. %. All mag- 
netic measurements were done using Quantum Design 
SQUID magnetometer (QDMPMS). The Curie tempera- 
tures (Tc) were determined from the rate of fall of mag- 
netization as a function of temperature in applied fields of 
20 - 50 Oe. The compositions of the samples have been 
confirmed by energy dispersive spectrum (EDS) measure- 
ments using a scanning electron microscope (SEM). The 
details of measurements of other quantities will be dis- 
cussed along with the presentation of data in section III. 



III. RESULTS AND DISCUSSION 
A. Magnetic phase of the alloys 

The X-ray diffraction pattern reveals that all the alloys 
are of single FCC phase with no additional phase. The 
lattice parameter varies by 1.5-2 % from that of pure Ni 
and increases with the increase of Mo/W. The alloys of 
both the Mo and W series show paramagnetic to ferro- 
magnetic transition at T c , which varies widely depending 
on the amount of Fe present in the alloy. The magnetiza- 
tion of the alloys with low Fe content are measured in ~ 
20 Oe applied field and with low temperature attachment 
of QDMPMS and that for the alloys with higher concen- 
tration of Fe were measured in ~ 50 Oe applied field with 
high temperature attachment of QDMPMS. The T c for 
all the compositions are already reported by Banerjcc ct 
al4£ and are given in Table I. It is clear that increas- 
ing Mo and W content and accompanying decreasing Fe 
content in the alloys lowers the T c rapidly. 

In the Mo series, the sample Mol3.5 shows a sec- 
ond phase transition apart from the paramagnetic- 
ferromagnetic transition found in all the alloys of the se- 
ries. No such additional phase has been found in the W 
series which we anticipate is due to lack of high W sam- 
ples around the critical concentration. We have found 
that the T c extrapolates to K for -14.2 at. % Mo 
whereas the Slater-Pauling curved gives Tc ~ K for 18 



% Cr in Ni. The split-band model of Berger— , however, 
gives the same critical concentration of — 12 % for Cr, 
Mo, and W since they fall in the same group in the pe- 
riodic table. Detailed magnetic measurements, both dc 
and ac, have been done on Mol3.5 to establish that the 
second phase is indeed a mixed ferro-spin-glass phase. 
The M(H) measurements show a ferromagnetic nature 
with a feeble non-saturating character at 5 K but grossly 
it is ferromagnetic. The peak in ac susceptibility shifts to 
higher temperatures with increasing frequency and also 
the field cooled (FC) and zero-field cooled (ZFC) curves 
bifurcate until an applied field ~ 100 Oe. These charac- 
teristic features resemble those of spin glasses. All the 
above will be presented in detail elsewhere^. 

The alloys of both the series show very thin M-H 
loops just like their parent permalloy. However, the per- 
meability found is not comparable to those of permal- 
loys. Magnetic annealing of the alloys could give higher 
permeability^. The saturation magnetization falls from 
0.85 hb to 0.08 [Ib on increment of Mo in the alloys^ 
while it decreases from 0.85 [Ib to 0.11 \lb for the W 
series^. The ac-susceptibility measurements show the 
same Tc as found from dc measurements. 



B. Low-temperature high-field magnetization of the 
alloys 

In conventional ferromagnets, the temperature depen- 
dence of the spontaneous magnetization M S {T) far be- 
low Tc is dominated by long-wavelength spin-wave ex- 
citations. The excitation energy E(k) of spin waves in 
the limit of small wave vectors (fc<Ca _1 , a is the lattice 
spacing) is given by 

E(k) = ftw(fc) = gfi B H iat + Dk 2 +Ek 4 + ..., (1) 

where the first term is an energy gap due to the pres- 
ence of an effective internal field i?i n t, arising from the 
applied field, the anisotropy field, and the spin-wave de- 
magnetizing field. D is the spin-wave stiffness constant, 
and E is a proportionality constant. Even in disordered 
ferromagnets, there is ample experimental evidence that 
long- wavelength spin- wave modes are a useful way to rep- 
resent the low-energy magnetic excitations^ 

In the low-temperature limit, according to the Heisen- 
berg model, the change in the spontaneous magnetization 
due to the excitation of spin waves can be written asJ^ 

M(T) = A/(0)[1-BzQ,^)t 3 / 2 

- Cz^Y^ + ...], (2) 

where M(0) is the magnetization at K, T s =g/i B if int /fc B 
is the gap temperature, and z ( |, -^-j and z ( |, -t^-J are 
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the correction terms which reduce to unity if the effective 
internal magnetic field vanishes. The two temperature 
terms above come from the harmonic (k 2 ) and anhar- 
monic (/c 4 ) terms in the spin- wave dispersion relation [Eq. 
(fT])], For simplicity, disregarding anharmonicity and the 
gap corrections, Eq. ([2]) reduces to the so-called Blochs 
T 3 / 2 law: 



M(T) = M(0)(1 - BT 3/2 ). 



(3) 



In conventional spin- wave theory the spin wave parame- 
ter B and the spin wave stiffness coefficient D are related 
through 



fee ( 2M2g^ B 
4tt ^ M(0)B 



2/3 



(4) 



The thermal demagnetization process of ferromagnetic 
metals at low temperatures (T<CTc) can be explained by 
both localized 1 ^ and itinerant^ models. In the localized 
model, with static interactions between spins associated 
with localized electrons on atomic sites, spin waves cor- 
respond to coherent superposition of local spin devia- 
tions in the ferromagnetic state, and their equilibrium 
number density at finite temperature according to the 
Bose-Einstein distribution function yields the magneti- 
zation reduction according to Eq. ([2]). In the itiner- 
ant model, electrons move in the average field of other 
electrons/ions, and effective inter-site spin couplings are 
generated by exchange of the particlc-holc propagator, 
strong correlation effects in which have been investigated 
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FIG. 1: (Color Online) Magnetisation vs. Temperature plot 
for four samples are given, among them three samples fit bet- 
ter using Eq. ((2j) including anharmonic term another one 
shown fitted to Blochs T 3/2 law. 



recently using a systematic non-pcrturbative expansion 
schemed The presence of thermally excited long wave- 
length spin waves in the ferromagnetic state at low tem- 
perature allows virtual excitation of majority-spin elec- 
trons to the minority-spin band due to electron-magnon 
coupling, resulting in spectral- weight transfer and conse- 
quent reduction of the magnetization (n-j-) — (^4.)- 

The low-temperature magnetization was measured 
from 2 K to 0.1 T c or till 15 K at ~ 6000 Oe for all the 
samples which is above their saturation fields of 1000- 
2000 Oe, although sample Mol3.5 did not actually satu- 
rate even at 6000 Oe due to its low-temperature mixed 
ferro-spin-glass phase. Figure 1 is a plot of M vs. T for 
T <C Tc. The solid lines are the best-fitted curves as in- 
dicated in the legends of the figures. The gap correction 
[Eq. (f2])] was tried but it was found to be negligible since 
the gap temperature came out to be less than 1 K. 

From Fig. 1 as well as Table I we find that the fits 
for some of the 12 samples that we have studied improve 
significantly if we include the T 5 / 2 term in Eq. (T2J) . Not 
only the values of \ 2 arc smaller by a factor of 2-3 (~8 
for W7.2), the deviation of the best-fitted curve from 
the data vs. T is random as against systematic when 
we consider only the T 3 ' 2 term. For the other samples 
the improvements are insignificant. We must note that 
for higher Tc samples (low Mo and W), AM/M is only 
1% and only very high resolution SQUID measurements 
are able to isolate the anharmonic term in the magnon 
dispersion relation. 

The first significant feature emerging from the present 
investigation is the order of magnitude of the spin-wave 
parameter B. These values are found to be strongly en- 
hanced with increasing Fe dilution in both the series of 
alloys, and are about 10 to 100 times larger than those 
found for conventional 3d fcrromagnets (for bulk Fe it is 
3.4 x 10~ 6 K~ 3 / 2 ). This sharp increase in the values of B 
with increasing Fe dilution reflects an enhancement in the 
density of low-energy magnetic excitations due to weak- 
ening of the ferromagnetic couplings between the bulk Fe 
spins forming the percolating ferromagnetic matrix. The 
spin-wave stiffness constants arc also correspondingly re- 
duced, as seen in Table I. 



C. Intermediate-temperature behaviour 



In order to allow a qualitative comparison between the 
magnetization behaviour of different ferromagnetic sys- 
tems over a much broader temperature scale extending 
nearly up to Tc, it is convenient to use the normalized 
coefficient B 3 / 2 defined through the relation i 18 ' 19 



M s (0)-M s (T) 
M s (0) 



B3/2 I ~ 



T 



3/2 



(5) 
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TABLE I: Magnetization vs. temperature data are fitted to Eq. ([2]) and the coefficients B, C along with normalized x 2 , R 2 , 
Tc, and D are tabulated below. 



Sample 


M(0) 


B 


C 


B3/2 


Tc 


normalized x 


t>2 


D 


Composition 


(emu) 


(10 -4 K _3/2 ) 


(10 -6 K _5/2 ) 




(K) 


(10 -7 ) 




meVi 2 


Nisi.oFei6.7Mo2.3 


3.33323±0.00005 


0.1189±0.0004 




0.23 


720 


0.099 


0.99843 


219 




3.33286±0.00004 


0.094±0.001 


0.023±0.001 






0.036 


0.99941 




Ni8O.4Fe12.gM06. 7 


1.4606±0.0002 


0.294±0.005 


- 


0.32 


495 


2.91 


0.97166 


144 




1.4601±0.0003 


0.24±0.03 


0.1±0.06 






0.036 


0.99941 




Ni 8 3.4Fei .7Mos.9 


0.79083±0.00006 


0.345±0.003 


- 


0.35 


470 


2.08 


0.99064 


155 




0.79038±0.00007 


0.24±0.01 


0.19±0.02 






1.28 


0.99439 




Nis3.5Fe7.6M0g. 9 


3.13870±0.00002 


0.864±0.007 


- 


0.49 


320 


3.31 


0.99981 


110 


Ni 8 3.iFe 6 .oMoi .9 


1.20957±0.00005 


3.68±0.01 


- 


0.90 


182 


0.75 


0.99916 


45 


Ni83.2Fe3. 3 Moi3.5 


0.2770±0.0001 


42.7±0.1 




1.29 


45 


63.9 


0.99871 


20 




0.2755±0.0001 


36.2±0.5 


0.44±0.03 






27.7 


0.99944 




Ni78.9Fei8.1W3.c1 


4.0269±0.0001 


0.0982±0.0007 




0.21 


775 


0.52 


0.99332 


255 


Ni79.4Fei4.iW 6 .5 


1.25119±0.00002 


0.302±0.001 




0.37 


530 


0.15 


0.99945 


133 




1.25098±0.00002 


0.275±0.002 


0.040±0.002 






0.07 


0.99976 




Ni8i.oFeii. 8 W 7 .2 


1.10312±0.00002 


0.359±0.001 




0.42 


515 


0.31 


0.99921 


137 




1.10280±0.00001 


0.311±0.002 


0.070±0.002 






0.04 


0.9999 




Ni 8 3.5Fe7.6W 8 .9 


0.56687±0.00001 


0.796±0.001 




0.41 


300 


0.47 


0.99956 


120 


Nis2.6Fe6.9W10. 5 


1.66087±0.00003 


2.906±0.004 




0.77 


191 


0.22 


0.99963 


53 


Nis6.eFe3.1W10. 3 


0.54141±0.00001 


19.6±0.1 




0.8 


55 


18.08 


0.99825 


28 



in terms of the low-temperature magnetization (Eq. [3j, 

3 /2 

which yields B^i 2 =BT c ' . Significantly, this reduced co- 
efficient B 3 / 2 provides an effective measure of the overall 
shape of the magnetization decay. For crystalline ferro- 
magnets, where the magnetization falls sharply near Tc, 
one obtains B3/2 ~ 0.2, irrespective of the Curie temper- 
ature. On the other hand, for reference, if the magnetiza- 
tion were to fall off as M(0)(1 - BT 3 / 2 ) nearly upto T c , 
then clearly -B3/2 ~ 1. Table Q] shows that with increasing 
Fe dilution, the -B3/2 values for both series of alloys sys- 
tematically increase from ~ 0.2 to ~ 1, indicating that 
the magnetization decay near Tc becomes significantly 
slower in the Fe-poor alloys. This is the second significant 
result which emerges from the intermediate-temperature 
regime. 

Both these distinctive features of macroscopic magne- 
tization behaviour observed in Ni-Fe-Mo and Ni-Fe-W 
alloys — strong enhancement of B with increasing Fe 
dilution and -B3/2 approaching unity — have also been 
observed in squid magnetization studies of diluted mag- 
netic semiconductors Gai-^Mn^Asji where they were as- 
cribed to two distinct energy scales involved in the local 
spin dynamics associated with the formation of spin clus- 
ters. A brief review of the theoretical analysis of finite- 
temperature spin dynamics within a minimal model for 
diluted magnets will be helpful in understanding the ob- 
served magnetization behaviour. 

Theoretical investigations in diluted magnets do 
yield competing ferromagnetic and antiferromagnetic 
interactions.— Within a minimal model for diluted mag- 
nets involving spin-S localized impurity spins and host 
band fermions (carriers), magnetic interactions between 



two impurity spins at lattice sites i and j were calcu- 
lated from Jij=J 2 (2S)4>ij in terms of the particle-hole 
propagator <pij evaluated in the ferromagnetic state. For 
the same impurity separation, the calculated magnetic 
couplings were found to exhibit a broad distribution, im- 
plying that the coupling between two impurity spins is 
not simply a function of their separation, but actually 
depends on the whole disorder configuration, suggesting 
shades of a complex system. 

The broad distribution in the calculated ferromagnetic 
spin couplings was ascribed to the formation of small 
impurity spin clusters due to positional magnetic disor- 
der. The cluster spin couplings were found to be strongly 
enhanced due to preferential accumulation of carriers in 
impurity spin clusters, whereas the consequent depletion 
of carriers from the bulk resulted in weakened bulk spin 
couplings between the bulk spins forming the ferromag- 
netic matrix^— 

The effects of such a broad distribution of ferromag- 
netic spin couplings on finite temperature spin dynamics 
was recently investigated in diluted magnets using a lo- 
cally self-consistent magnon renormalization scheme^r— 
in analogy with renormalized spin- wave theory in ordered 
ferromagnets^ The local magnetization (Sf) of each in- 
dividual impurity spin at site i in a quantum spin-S disor- 
dered fcrromagnet was obtained by self-consistcntly solv- 
ing the three coupled equations: 

, z] (S - + $,) 2A ' +1 + (S + 1 + $ t )$f +1 

{ i ' (1 + $,)2S+1 _ ^S+l 

(6) 

where the local site-dependent boson occupation num- 
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bers: 

d>.=V^- (7 ) 

z 

in terms of the renormalized magnon energy eigenvalues 
oji and eigcnfunctions 4>i obtained from the renormalized 
magnon Hamiltonian: 

Uij = y/2{sf) (J 2 [x%) y/2(sf} ■ (8) 

The locally self-consistent magnetization (5f } thus ob- 
tained clearly showed rapid thermal demagnetization and 
nearly paramagnetic behaviour of the weakly coupled 
bulk (FM matrix) spins, whereas the strongly coupled 
cluster spins resist thermal demagnetization and thus 
prolong the magnetic order near Tq- The overall pic- 
ture from the averaged magnetization was that while the 
ferromagnetic Xc was suppressed by positional-disorder- 
induced spin cluster formation, the magnetization decay 
was distinctly stretched near Tq. 

These distinctive features of microscopic spin dynam- 
ics behaviour in a generic diluted magnet were shown 
to quantitatively affect the two readily accessible char- 
acteristics B and B 3 / 2 of macroscopic magnetization be- 
haviour. While the calculated spin- wave parameter B 
was found to be sharply enhanced with increasing dilu- 
tion [Fig. 2 (upper panel)] due to weakened bulk spin 
couplings and softening of low-energy spin excitations, 
the stretching of magnetic order near Tq due to strongly 
coupled cluster spins resulted in enhanced B 3 / 2 ~l, in- 
dicating slower thermal demagnetization approximately 
as T 3 / 2 (dashed line) over a much broader temperature 
range, as shown for five different disorder configurations 
(lower panel). In sharp contrast, the ordered ferromagnet 
yielded, for the same dilution and carrier concentration, 
a conventional magnetization decay with -B3/2 ~ 0.2 and 
T 3 / 2 behaviour (dashed line) only in the low-temperature 
regime. Here the given notation refers to impurity con- 
centration (x), carrier concentration (p), bandwidth (W), 
and host-impurity energy offset (e^). 

As shown in Table HI a very similar behaviour is ob- 
tained for the two magnetization (spin-wave) coefficients 
B and -B3/2 from our finite-temperature magnetization 
investigation of Ni-Fe-Mo and Ni-Fe-W alloys. While the 
B values are sharply enhanced with increasing Fe dilution 
(the spin stiffness is correspondingly sharply suppressed), 
the -B3/2 values approach I with increasing Fe dilution in 
both alloys, as compared to about 0.2 for ordered ferro- 
magnets. Similar enhancements were reported recently 
from SQUID magnetization measurements on the diluted 
magnetic semiconductors Gai_ x Mn x As^ and in earlier 
studies on amorphous ferromagnetic alloys in compari- 
son with crystalline ferromagnetsi 14 ' 15 
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FIG. 2: (Color Online) Rapid enhancement of the calcu- 
lated spin- wave parameter B with dilution (upper panel), and 
(lower panel) nearly T 3 ^ 2 fall-off of site-averaged magnetiza- 
tion and stretching of magnetic order near Tc (from ref. [7]). 



D. Magnetic phase transition and the associated 
critical exponents 

In the present work, we have also studied the criti- 
cal exponents associated with the magnetic phase transi- 
tion of the alloys containing more than 8 at.% of Mo/W. 
The specific reason to study the critical exponents and 
the amplitudes is to see the consistency of our findings 
with respect to the above magnetization data analysis. 
In the last section we have seen evidence of spin clusters 
with relatively strongly coupled spins in the alloys with 
dearth of Fe content. How does this spin clustering af- 
fect the critical exponents and critical amplitudes? As 
a first step, ac susceptibility has been measured in all 
alloys with a temperature increment of 0.2 K or less so 
that their Tc could be found accurately within ± 0.1 K. 
The magnetization was measured in fields ranging from 
to 20 kOe at various temperatures around 2% of Tc on 
either side. 

The critical exponent (3 associated with M s , the spon- 
taneous magnetization, is given by M s =B\e\@ for T<Tq, 
where e—(T — Tc)/Tq and B is the corresponding crit- 
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ical amplitude. The critical exponent 7 is related to 
Xo, the zero- field dc susceptibility, through the equation 
Xo 1= r _1 |£| 7 , where T _1 is the critical amplitude. The 
magnetic field dependence of M at Tc gives us the third 
critical exponent 8 through the relation AI=DH X / S , D 
being the critical amplitude. All these three exponents 
follow a static scaling relationship, 8=1 + 7//?. In order 
to find these exponents one has to know the transition 
temperatures very accurately. First the isothermal mag- 
netization M was plotted as in a simple Arrott plot (M 2 
vs. H/M) which uses the mean-field exponents. How- 
ever, when they did not give linear isotherms, the modi- 
fied Arrott- Noakes (AN) plot was tried, i.e., M 1 ^ versus 
{H/M) 1 /" 1 . 

To plot the AN isotherms one has to know the values of 
/3 and 7 beforehand. Instead of guessing the initial values 
we have calculated j3 and 7 using a simple FORTRAN 
code where we applied the method of least squares fit- 
ting. We have calculated the average chi-square (x 2 ) for 
all possible combinations of /3 and 7, where /3 was var- 
ied between 0.2 and 0.6 in steps of 0.001 and 7 from 
0.9 to 1.5 in similar steps. Then the combination of 
/3 and 7, for which the average x 2 is a minimum and 
the slopes of all the straight lines are almost equal, i.e., 
the set of isotherms are truly parallel to each other, are 
taken. Amazingly, the extrapolated isotherm at Tq in- 
deed passes through the origin making us quite confident 
of the plots. Then using those values of /? and 7 we have 
plotted in Fig. [3] the modified Arrott- Noakes (AN) plot, 
a set of isotherms for 10 % Mo alloy, the y-intercepts 
of these straight lines give M s and the x- intercepts Xo ■ 
Figure [3] (upper panel), a typical AN plot, also gives the 
values of /? and 7 for the best-fitted isotherms. Although 
7 is close to the mean- field value of 1, (3 is much lower 
than the mean-field value of 0.5. Deviations of the values 
of j3 and 7 from the mean-field values have been found 
both in many magnetic glasses as well as in crystalline 
fcrromagnets. 

There is another way of calculating j3 and 7, that is 
using the Kouvel-Fisher method. The set of equations 
used in this method are: 



2.0x10" 



Y{T) 



/d\nM s 
V dT 



{T - T C ) (9) 



and 



X{T) = 



MnXo 1 
^ dT 



= (~)(T-T C ). (10) 



These equations arc valid only in the critical region where 
temperatures are very near Tc. Here we have used the 
M s and XcT* values from the intercepts of the AN plots. 
In the critical region Y{T) vs. T and X{T) vs. T are both 
straight lines with slopes of (1//3) and (1/7), respectively. 
In this method a priori knowledge of Tq is not needed and 
if the results are consistent, then Y{T) and X{T) will 
intersect the T axis at the same point. Figure [3] (lower 
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FIG. 3: (Color Online) The modified Arrott-Noakes (AN) plot 
and X (T) and Y (T) vs. temperature plots for Mol0.9 alloy 
with values of ft and 7 for which the best-fitted isotherms are 
obtained. 



panel) also shows the typical X{T) and Y{T) plots for 

10 % Mo alloy. Both of them intersect the T-axis at 
T=T C =180.5 K. 

We have used both the methods for all our samples to 
check the consistency of the results. Except for Mo 13. 5 
alloy, all the other samples gave quite reasonable val- 
ues of the critical exponents. However, it is to be noted 
that unlike other alloys, Mo 13. 5 has a reentrant spin-glass 
phase (RSG) phase below T 9 =10 K (T c is ~ 44 K) which 
might possibly affect the critical behavior since the tran- 
sition may not be a pure ferro-para one although we used 
M (T, H) data only from T=43-45 K, far away from 10 K. 
Nevertheless, our calculated values of the critical expo- 
nents from the two methods did not match at all. More 
importantly, the value of (3 was above and that of 7 was 
below the mean-field values. This is quite unphysical as 
the mean-field case is the limiting one, since it is a rather 
crude theory. Presence of more than one phase due to 
improper homogenization and/or final annealing might 
cause such deviations in the critical exponents. Table 

11 below gives the values of Tc, /?, 7, and 8. The criti- 
cal exponent values suggest that Mo8.9 and W8.9 almost 
follow the mean-field model whereas the others, except 
Mol3.5 are close to 3D-Ising rather than that of pure 
nickel, which follows the 3D-Heisenberg model. 

We have also calculated the critical amplitudes since 
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TABLE II: Alloy compositions, values of Tc and critical ex- 
ponents, obtained both experimentally and from KF analysis 
along with those for pure Ni and those from existing theories. 



Sample 


T C 


P 




7 




5 


1 


Composition 


Expt. 


KF 


Expt. 


KF 


Expt. 


KF 


Expt. 


Nig3. 5 Fe7.6Mo 8 .9 


316.0 


316.4 


0.48 


0.42 


1.21 


1.31 


3.54 


3.50 


Nis3.iFe6.oMoio.9 


180.5 


180.4 


0.275 


0.31 


1.03 


1.05 


4.68 


4.74 


Ni 8 3.2Fe3. 3 Moi3.5 


44.55 


46.2 


0.59 


0.54 


0.92 


0.88 


2.22 


1.64 


Ni83.5Fe7.6Ws. 9 


305.0 


300.0 


0.43 


0.45 


1.30 


1.37 


4.22 


4.05 


Nig2.6Fe6.9W10. 5 


191.0 


190.0 


0.34 


0.39 


1.23 


1.19 


4.67 


4.64 


Nig6.6Fe3.1W10. 3 


58.15 


58.15 


0.33 


0.35 


1.304 


1.29 


4.93 


4.54 


Ni a 


627.4 




0.378 




1.34 




4.58 


4.54 


Mean-field 






0.50 




1.00 




3.00 




3D-Ising 






0.312 




1.25 




5.00 




3D-Heisenberg 






0.378 




1.405 




4.76 





Reference 21. 



4 Ni. 



Fe W 

.6 3.1 10.3 



y 

y = 1.3 

= 4.63 kOe g/emu 



NL 



4.9 4.2 

In e 
, Fe ,o M %., 



■3.5 



: 0.289 

= 19.95 emu/g 



4.9 4.2 
In £ 



Ni »,e Fe e 



= 4.67 
= 5.13 



8.4 
InH 



NL.Fe. W 



86.6 3.1 10.3 



10 

10 : 



H/|s|P 5 



they are important for a complete knowledge of the 
critical behavior near ferro-para transition. The mag- 
netization above and below Tc satisfies a single scal- 
ing equation given by m=f±(h) where m=\e\~ @ M (e , H) 
and h=\e\~P s H are called scaled magnetization and the 
scaled magnetic field. The above relations show that m 
as a function of h falls on two different universal curves 
/_0) for T<T C and f+(h) for T>T C . If the values of 
the critical exponents found here are correct, then all the 
data will fall on two distinct curves confirming their cor- 
rect choice. The critical amplitudes Bo=mo, T~ 1 =ho/nio 
and Dq are obtained from the intercepts of the ln-ln plot 
of M a vs. |e|, Xo 1 vs - l £ l an< ^ M vs. H(e = 0), re- 
spectively. Figure S] shows, respectively lnx^~ vs. lne 
plot for W10.3, In M s vs. InH for W10.5, lnM s vs. lne 
for Mol0.9, and M/\ef vs. H/\e\ 05 for W10.3 including 
the calculated values of the corresponding critical am- 
plitudes. Figure [H (lower right panel) shows clearly the 
correctness of our procedure. Table IIIII gives the criti- 
cal amplitudes of some of the samples. We observe that 
with the increasing Fe dilution, the critical amplitudes 
r _1 , Bo, and Do decrease systematically and all three 
have values lower than those for Ni. This is the third 
observable consequence of strongly coupled cluster spins, 
indicating significantly reduced participation by the bulk 
spins forming the ferromagnetic matrix in the critical be- 



TABLE III: The critical amplitudes of the alloys and those of 
Ni for comparison. 



Sample 


r- 1 


Bo 


Do 


composition 


(kOe-g/emu) 


(emu/g 




Nig3.5Fe7.6M0g. 9 


8.4 


21.4 


7.36 


Ni 8 3.iFe 6 .oMoio.9 


4.2 


20.0 


5.14 


Nig3. 2 Fe3. 3 Moi3.5 


1.7 


2.5 


1.9 


Nig3. 5 Fe7.6Wg.9 


8.7 


27.2 


7.27 


Nig2.6Fe6.9W10. 5 


6.1 


25.9 


5.13 


Nig6.eFe3.1W10. 3 


4.6 


8.8 


4.79 


Ni a 


19 


83 


30 


a Reference 21. 



FIG. 4: (Color Online) . Xo* vs - m£ P lot for W10.3 alloy, 
In M s vs. InH plot for W10.5 alloy, lnM s vs. lne plot for 
Mol0.9 alloy and M/\e\ p vs. H/\ef s plot for W10.3 alloy. 



haviour. 

The locally self-consistent spin dynamics analysis 
shows that as the temperature falls below Tc, the cluster 
spins rapidly get magnetized^— A similar behaviour is 
expected above Tc in an external magnetic field. Thus, 
the cluster spins contribute dominantly to the high mag- 
netic susceptibility of the alloys near Tc , and hence result 
in a suppression of Xo 1 an< ^ °f the critical amplitudes. 

Particularly for the Mol3.5 sample, the extremely low 
measured values of critical amplitudes along with the 
unusual critical exponent values (see Table II) indicates 
dominant finite-cluster contribution just below the per- 
colation threshold, with no diverging spin correlation 
length and no true critical behaviour. This is consistent 
with the emergence of spin-glass behaviour at this com- 
position due to frustrated cluster spins locked in frozen 
non-collinear orientations which also results in smaller 
spontaneous magnetization along the z-direction. 



IV. CONCLUSION 

In conclusion, we find that all three macroscopic mag- 
netization characteristics — the spin- wave parameter B, 
the reduced coefficient -B3/2, and the critical amplitudes 
T _1 — corresponding to the low, intermediate, and the 
critical temperature regimes respectively, yield distinc- 
tive spin dynamics signatures associated with strongly 
coupled cluster spins in these alloys. 

The dilution behaviour of the two magnetization co- 
efficients extracted from the macroscopic magnetization 
behaviour of finite-temperature spin dynamics in the two 
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alloy systems is indicative of a broad distribution of mag- 
netic spin interactions between the magnetically active 
Fe atoms. The sharp enhancement in the measured spin- 
wave parameter B accompanying the spin stiffness reduc- 
tion with Fe dilution clearly indicates weakened bulk spin 
couplings and softening of low-energy spin-wave modes. 
Furthermore, the magnitude of the reduced coefficient 
.B3/2 rapidly approaches 1 with increasing Fe dilution, 
and the slower magnetization decay with temperature in- 
dicates presence of strongly coupled cluster spins which 
resist thermal demagnetization and stretch the magnetic 
order near Tq. 

As also observed in diluted magnetic semiconductors 
and metallic glasses, this sharp enhancement with dilu- 
tion is thus suggestive of the two magnetization coeffi- 
cients as universal macroscopic indicators of spin cluster- 
ing and disorder-induced distribution of magnetic inter- 
actions in disordered ferromagnets. The reduction of the 
measured critical amplitudes with Fe dilution due to sig- 
nificantly reduced participation by the bulk spins form- 
ing the ferromagnetic matrix in the critical behaviour is 



consistent with this picture. Similar suppression of criti- 
cal amplitudes was observed in Fe x Niso- x Pi4B6 alloys^ 
and was ascribed to the growth of spin clusters as T ap- 
proaches Tq and the consequent reduced participation of 
remaining bulk (FM matrix) spins in the FM-PM phase 
transition. 

Except for the Mol3.5 alloy, which has a spin-glass 
phase below T g =10 K, all the other samples gave 
quite reasonable values of the critical exponents. Non- 
universal values of critical exponents with j3 « 0.55 have 
also been found in recent experimental studies^ of the 
ferromagnetic transition in spin-glass re-entrant metallic 
alloys such as Auo.8iFeo.19. 
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